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ABSTRACT 

Context. There is mounting evidence that a significant fraction of Black Holes (BHs) today live in late-type galaxies, including bulge- 
less galaxies and those hosting pseudobulges, and are significantly undermassive with respect to the scaling relations followed by 
their counterpart BHs in classical bulges of similar stellar (or even bulge) mass. 

Aims. Here we discuss the predictions of two state-of-the-art hierarchical galaxy formation models in which BHs grow via mergers 
and, in one, also via disk instability. Our aim is to understand if the wealth of new data on local BH demography is consistent with 
standard models. 

Methods. We follow the merger trees of representative subsamples of BHs and compute the fractional contributions of different 
processes to the final BH mass. 

Results. We show that the model in which BHs always closely follow the growth of their host bulges, also during late disk instabilities 
(i.e., bars), produces too narrow a distribution of BHs at fixed stellar mass to account for the numerous low-mass BHs now detected 
in later-type galaxies. Models with a looser connection between BH growth and bar instability instead predict the existence of a larger 
number of undermassive BHs, in better agreement with the observations. 

Conclusions. The scatter in the updated local BH-bulge mass relation (with no restriction on galaxy type) appears to be quite large 
when including later-type systems, but it can still be managed to be reproduced within current hierarchical models. However, the 
fuelling of BHs during the late bar-instability mode needs to be better quantified/improved to properly fit the data. We conclude 
discussing how the possibly large number of BHs in later type galaxies demands for an in-depth revision of the local BH mass 
function and its modelling. 

Key words. Galaxies: bulges - Galaxies: evolution - Galaxies: nuclei - Galaxies: statistics - Galaxies: structure - Cosmology: theory 



1. Introduction 

There is evidence in the local Universe that most bulge- 
dominated galaxies seem to obey rather tight correlations 
between the central black holes' (BHs) masses and the 
. . . bulge masses and velocity dispersions of the ir hosts (e.g., 
> ■ |Ma gorria n et al . 1998: lFerrarese & M erritt 20Q3). 

• i-H ■ 

Such tight correlations have been physically interpreted as 
5_i " byproducts of some BH self-regulation mechanism - active 
C3 , galactic nuclei (AGN) feedback - operating during the high- 
redshift, gas-rich phases of galaxy build-up. At early times in 
fact, large amounts of gas inflows can be gen erated due to, 
e.g., merger events o r internal instabilities (e.g., lBarnes|[l992t 
Bournaud et al.ll20TTh . triggering star formation and fuelling the 
central BH that can energetically and kinematically back-react 
on the gas reservoirs via AGN winds. Such interactions can 
possibly contribute to shut-down starformation and to the for- 
mation of the progenitors of today 's earl y-type galaxies (e.g. , 
Sanders et al] 1 1988t iGranato et all l2006t IShankaret ail [2008; 
Lapi et al.ll201 lb . 



that BH-host correlations indeed appear as tight as previously 
claimed for all typ es of galaxies (e.g., iGultekin et al.1 2009; 



Mathu r et alj 



Mounting empirical evidence suggests, however, that other 
additional tracks may have significantly contributed to BH 
formation and growth. Different observations now question 



i Graham & Lil2009l). Significa nt outliers exi st (e.g ., 
GOO It iGrupe & Mathurfooofl IWatson etafll2007l) . BHs are cur- 
rently observed also in the centres of bulge-less galaxies (e.g ., 
iGhosh et alj2008HSatvapal et al.l2009HArava salvo et al.l2 012), 
and the AGN-starformati on correlations disappear for all b ut 
highest luminosity AGNs (IShao et al . 201ft iGrier etaTEoilh . 

Perhaps one of the most pressing reason to argue for supple- 
mentary routes of BH e volution is the prevalence of pseudob- 
ulges hosting BHs (eg.,lHdl2008LlGadotti & KauffmannlF2 009 ; 
I Kormendv et al.ll201 UlMathur et al.ll201 lUOrban de Xivrv et al.1 
1201 lb . "Pseudobulges" are structurally and dynamically dif- 
ferent from "classical" bulges and ellipticals, which can also 
have very low specific angular momentum (e.g., Emsellem et al.l 
2011). They appear to have more "disky" profiles, they depart 
from the fundamental plane of elliptical gala xies, and are fainter 
at a fixed effective radius (e.g.,|Gadotti 200§). 

At variance with elliptical galaxies and classical bulges, 
pseudobulges have been suggested to h ave formed through sec- 
ular p rocesses such as disk instabilities (IKormendv & Kennicuttl 
l2004i) . What we refer to here is bar instability, physically distinct 
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Fig. 1. Le ft panel: The dark blue circles show the distribution of black holes for a sample of 1000 black holes randomly selected 
from the iBower etal ] (l2006h catalogue. The cyan squares are a random sample of 100 black holes with their hosts' bulges mainly 
grown via secular processes ("pseudobulges"). The long-dashed yellow and red black lines are the median Mbh-Mb u i ge relations 
competing to the for mer and the latter samples, respectively; the grey stripe indicates the local empirical Mbh-Mb u i ge relation for 
classical bulges from ISani et al.l (1201 ll) . Middle panel: same pattern as the left panel but with both samples extracted from the 
iGuo et all (1201 ll) catalog, with the long-dashed yellow and red solid lines being the respective median Mbh-Mbuige relations for the 
two samples. Right panel: comparison between predicted the distributions of pseudobulges of IBower et al.l d2006[) and Guo et al. 
(201 1; bl ue and red lines, respect ively), with thei r 1-cr uncertainties (dashed lines) com pared to data (Table 1 ) from | Hul (I2009L filled 

1\Q~ 



squares), Kor mendv et al.l (1201 ll filled circles). iMathur et ail ( 201 ll ope n circles), and lOrban de Xivrv et alj d201 ll open squares). 
Also shown for completeness data on pseudobulges bv lSani et alh 201 ll open triangles). The average error in BH masses for active 
and inactive galaxies is of the order of ~ 0.3 dex. Of similar magnitude is the propagated error on stellar masses at fixed galaxy 
colour. 



to the clu mpy instabilities that characterize the h igh-z, gas-rich 
disks re.g.. lbaddi et al.ll2010tlGenzel et al.ll20flh . 

Indeed, [Shankar et al.l d201 lal) have recently shown that 
bulges m a inly f ormed from bar-like disk instability in the 
IGuo et ail (1201 ll) model have much lower bulge-to-total (B/T) 
ratios and smaller sizes t han classical bulges at fixed total stel- 
lar mass. iHopkins et a ll d2010h have argued that the formation 
of pseudobulges dep ends on how gas-rich the merger is, and 
Font anot et alj d201 lb claim that it depend s on the mass of the 
dar k matter halo (see also discussions in iFanidakis et alJl201ll 
and lYates et al.ll2"oTlt) . In particular, as detailed below, BHs re- 
siding in pseudobulges appear to be relatively undermassive 
compared to those in cl assical bu lges of similar luminosity or 
velocity dispersion (e.g., Hu 2008). 

Despite the strong observational efforts carried out over the 
years to characterize the BH population in later-type galaxies, 
theoretical models have mostly focused their attention on BHs 
in classical bulges. Here we investigate the formation of late- 
type bulges and their nuclear BHs, with emphasis on the class 
of galaxies broadly referred to as pseudobulges, in the frame- 
work of state-of-the-art hierarchical galaxy formation models 
that grow bulges and/or their central BHs via mergers and disk 
instabilities. In particular, we will show that mergers can produce 
a significant scatter in the BH-galaxy mass scaling relations, in 
broad agreement with the observed properties of classical and 
pseudo bulges. 



2. Reference Models 

Th e results discussed h ere are based on the on line catalogues of 
the IBower etal] d2006l B06 hereafter) and the lGuo etalfd2011l 
Gil hereafter) semi-analytic models Q. 

Both models follow the evolution of galaxies and their cen- 
tral BHs along the hierarchical evolution of their host dark mat- 
ter haloes and subhaloes within the concordance ACDM cos- 



mology. To this purpose, they are both implemented on top of 
the same large, high -resolution cosmolog ical N-body simula- 
tion MILLENNIUM I dSpringel et al.l l2005h . making their com- 
parison even more meaningful. The MILLENNIUM I simulation 
follows the evolution of /V = 2160 3 dark matter particles of mass 
8.6 x 10 8 h~ l M Q , within a comoving box of size 500/i _I Mpc on 
a side, from z = 127 to the present, with cosmological parame- 
ters Q m = 0.25, O b = 0.045, h = 0.73, Q A = 0.75, n = 1, and 
erg = 0.9. The simulation can resolve haloes down to a mass of 
> 10 10 M Q /h, corresponding to BHs close to the seed mass used 
in the models (> 1O 4 M ), thus, any conclusion drawn for BHs 
more massive than ~ 10 5 M is free of any resolution effect. 

The primary channel to form bulges in both models is via 
mergers, during which bulges gain a fraction of the progeni- 
tor stars along with newly formed stars, both dependent on the 
strength and type of the merger. 

Both models also substantially grow bulges via bar-like disk 
instabilities. Secular evolution is believed to form bulges via a 
bar that is capable of transferring mass from the disk to the cen- 
tral regions. At each timestep for each galactic system a stability 
criterion of the type 



Vmax < V, 



ref > 



(1) 



http://gavo.mpa-garching.mpg.de/MyMillennium3 



is checked. Here V max is the velocity at the maximum of the ro- 
tation curve, while V re f °c VGMdi s k//?disk is proportional to the 
actual velocity of the disk via a factor close to unity. 

When Eq.[T]is satisfied the models assume that the instability 
drive s the formation of a bar from the disk into a central bulge 
(e.g.. ICombes et alJll990b ICole et al.ll2000l) . Despite the broad 
scheme of the disk instability is apparently similar in both mod- 
els, there are significant differences. Gil treat disk instabilities 
as secular processes that act in self-gravitating disks near the in- 
stability boundary by transferring a portion of the stellar mass 
<5M stal from the disk to the bulge to keep the disk marginally sta- 
ble with an exponential density profile. 

B06 instead assume that the entire mass of the disk is trans- 
ferred to the galaxy bulge when the disk goes unstable, with any 
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Fig. 2. Black hole mass distri butions for galaxi es with bul ge masses 10 < log Mb u i ge < 10.5 {left panel) and 10.5 < log Mb u i ge < 1 1 
(right panel) predicted by the lGuo et alj d201 lb and lBower et al.l ((2006) models (red and blue lines, respectively), and for classical 
and pseudo bulges (dotted and solid lines, respectively). The solid, black lines in both panels are the data collected in Table Q] 
The grey stripes are the black hole mass intervals (median plus scatter) inferred from the Sani et al. Mbh-Mb u i ge relation for the 
corresponding bulge masses. 



gas present assumed to undergo a starburst and a fraction of it 
to feed the BH. It is thus expected that disc instabilities play 
a more prominent role in the B06 model. B06 indeed find that 
in their model most of the BH growth occurs via disk instabili- 
ties rather than mergers. The disk instabilities considered by B06 
may also qualitatively share some features with the instabilities 
mentioned in Sect . Q] characterizing high-z, gas-rich disks (e.g., 
iGenzel et alJ201 ll) . in the sense that in the B06 a significant frac- 
tion of gas can be fuelled to the centre even during disk instabil- 
ities. Nevertheless, the B06 are still considered and modelled as 
ba r-like instabilities accompanied by a starburst (see discussion 
in ICole et al.l 12000). thus still quantitatively different from the 
clumpy ones discussed in Sect. [T] Galaxy evolution models still 
lack a proper treatment of bulge formation (and any consequent 
BH feeding) via clump migration in gas rich , turbulent system s 
(see further discussions in Sect.[5]and in, e.g jDekel et alj2 009). 

Both models postulate that at each merger the central BH ac- 
cretes mass proportionally to the available gas in the host galaxy 
(in the Gil model the accretion efficiency gets lower for smaller 
mass systems and for unequal mergers), as well as mass from 
other incoming BHs. 

The most noticeable difference relevant to the present work 
is that, at variance with the Gil model, BHs in the B06 model 
also grow in mass during disk instabilities, proportionally to 
the fraction of gas available and undergoing a starburst at the 
moment of the disk instability. That some of the gas can be 
channelled into the very central regions during disk instabilities 
echoes some of the earlier models such as the so-called "bars 
within bars" schem e, proposed as a pos sible mechanism for fu- 
eling AGNs (e.g JShlosman et alJ l989). For suitable conditions 
the gas inflow induced by the bar instability may produce a self- 
gravitating disk in a smaller region and a new small bar can de- 
velop, which induces further gas inflow. 

In both models BHs also slightly grow in mass via hot ac- 
cretion at late times (radio-mode feedback), but this is a minor 
channel of BH growth compared to cold gas accretion and BH- 
BH mergers. Anyway, including or neglecting such processes in 
the models is irrelevant to the present discussion. 



Both models have been fine-tuned in their input parameters 
to reproduce the local scaling relations between BH mass and 
host galaxy stellar mass. However, the intrinsic scatter in the lo- 
cal relations is a genuine prediction of the models. The scatter 
reflects the relative diversity in the growth histories of central 
BHs hosted by galaxies of similar mass, induced by the intrinsic 
dispersion in the merger histories, and the different underlying 
physical assumptions (mergers and/or disk instability) to grow 
BHs in the models. Our discussion below will start from the 
comparison of the predicted scatter against the wealth of data 
now available on the BH-galaxy scaling relations. 

3. Data 

In the following, we will compare model predictions with a 
collection of data taken from the Literature, listed in Table 1. 
These data are for active and inactive BHs the hosts of which 
have been mostly recognized as pseudobulges. The column 3 of 
Table 1 lists references for both BH and bulge masses. Where 
stellar masses were not directly available, we obtained them 
by converting K, r, or V mag nitudes using the mass-colour re- 
lations from iBell et a D J20031) with an average g — r — 0.5, 
typical of these galaxies. All t he data and the models have 
been converted to a common JChabrierl (120031) Initial Mass 
Function. In lKormendv et al] (1201 ll) BH masses have been mea- 
sured using dynamical me thods. In iMathur et alj d201 ll) and 
lOrban de Xivrv et aH (1201 ll) widths of the H/3 emission lines in 
AGN spectra have been used for BH mass estimates (see Mathur 
et al. 201 1 for details. Orban de Xivry et al. do not list the BH 
masses in their paper; these have been kindly supplied by the au- 
thors). Hu (2009) and Sani et al. (201 1) report BH masses from 
multiple sources (please see their papers for details). 

Different authors have used different techniques to estimate 
BH (and bulge) masses, and each technique has its own set of 
errors. The average error in BH masses for activ e and inactive 
galax ies is of the order of ~ 0.3 dex (see, e.g., iMathur et al.l 
2011). Of similar magnitude is the propa gated error on s tellar 
masses at fixed galaxy colour (see, e.g.. IBell et al.l l2003h . In 
the following, for each galaxy we will use and show the mul- 
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Fig. 3. Comparison among the median properties of random su bsamples of classi cal and pseudobulges from the iGuo et al 1 (1201 lh 
model (red, dot-dashed and blue, solid lines, respectively) and lBower et al.l (feOOol) model ( red, long-dashed and black, solid lines, 
respectively). Upper panels: bulge mass (left), bulge-to-total B/T ratio (middle), and black hole mass (right) as a function of cosmic 
time or redshift. Lower panels: subhalo virial host mass (left), star-formation rate (middle), and specific star-formation rate (right) as 
a function of redshift. Overall, pseudobulges appear to have prolonged star-formation rates, significantly lower B/T, and in models 
with only mergers as fuelling mechanism, lower mass black holes, in agreement with observational data. 



tiple available BH and relative bulge mass estimates evaluated 
by each group. As it will be shown and discussed below, while 
some of the scatter in the data is undoubtedly induced by differ- 
ent measurement techniques, despite the different selections and 
methodologies, these data share a similar portion of the M^h- 
Mbuiee plane, mainly localized below the classical Mbh-Mb u i ae 
relation. The latter is not a new result and is discussed in the 
papers listed above. Here we simply put these data together for 
co mparison with model s. 

IGreene et ail d2008l) and iJiang et al.l d2011l) have confirmed 
the results of the a bove mentioned wo rks, extending them to 
lower BH masses. IGreene et al.1 (1201 Ol) also report lower BH 
masses at fixed velocity dispersion for a (smaller) sample of lo- 
cal late-type galaxies with more accurate BH masses from maser 
kinematics. However, we do not report their results in Table 1 
as only one galaxy in their sample with secure bulge masses is 
recognized as a pseudobulge. Also the indirect arguments put 
forward by Gadotti & Kauffmann (2009) suggest pseudobulges 
have BHs less massive than their classical counterparts at fixed 
bulge mass or velocity dispersion. 



4. Results 

4.1. A narrow and a broad distribution of Black Holes 

From now on we will simply refer to bulges grown mainly via 
bar-like disk instabilities as "pseudobulges", as opposed to "clas- 
sical" bulges that have grown their mass mainly via mergers. 
Such a classification is rather arbitrary, though refined theoret- 
ical and observational work is now establishing that movement 
of gas and stars via bars to the centre could effectively contribute 
in building a young and kinematically cold stellar bulge compo- 



nent, that is, a disc-like/pseudo bulge (e.g.. iGadottil |20 1 H and 
references therein). M oreover, it was recently pointed out by 
IShankar et al.l d201 lal) that galaxies of similar stellar mass that 
underwent a different merger history appear structurally differ- 
ent in the local Universe. In particular, they emphasized that in 
the Gil model those galaxies that had their bulges mostly grown 
via disk instability end up having significantly lower bulge-to- 
total ratios B/T and lower h alf-mass radii at fixed stellar mass, 
in agreement with data (e.g^ Gadotti 2009). Here we extend the 
differences emphasized bv lShankar et al.l (120 
and other properties. 



12031 to BH masses 



The blue circles in the left panel of Fig. 1 are a random 
subsample of 1000 BHs extracted from the B06 catalog, with 
the yellow, dashed line marking the median relation. A rela- 
tively narrow Mbh-Mb u i ge relation is predicted in fair agreement 
with what empirically calibrated in the local Universe (the grey 
band is the Mbh-Mb u i ge fit for classical bulges by Sani et al. 
2011 with its intrinsic s catter; similar results were found by, 
e.g., lHaring & Rixl l2004. though the latter are mainly based on 
dynamical masses). The cyan solid squares represent 100 ran- 
domly extracted BHs hosted by pseudobulges, with the red, 
long-dashed line marking the median relation. These galaxies 
were initially selected from the catalogs as "potential" pseudob- 
ulges for having lower B/T (usually set to be B/T< 0.5) at fixed 
bulge stellar mass. By analyzing their merger trees, following 
Shan kar et al.l (1201 lal) we then selected and defined as pseudo 
(classical) bulges those galaxies that had grown more (less) than 
50% of their bulge mass via disk-instability. We repeated the ex- 
ercise of randomly extracting 100 galaxies with such character- 
istics several times finding similar results. Given that BHs even 
in disk-instability mode in the B06 model grow proportionally 
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to bulges, it is not unexpected that they will end up following 
similar scaling relations to the BHs in classical bulges. 

The middle panel of Fig. 1 shows the predictions of the Gil 
model. The predicted Mbh-Mb u i ge relation (blue circles), with the 
median shown with a yellow, long-dashed line, is still broadly 
consistent with local data, though the predicted scatter at fixed 
stellar mass is much larger in this model. The cyan squares are 
a subsample of 100 randomly selected BHs that have grown 
more than 50% of their bulge mass via bar-like disk instabil- 
ity. Galaxies of the latter type were again initially selected as 
having lower B/T at fixed bulge mass. It is clear that extract- 
ing merger histories for the millions of BHs in the online cat- 
alogues becomes prohibitive. Nevertheless, for the Gil model 
we compared our randomly sel ected pseudobulge sa mple with 
one directly extracted from the Shan kar et al.l d201 lab catalog. 
The latter is a rendition of the Gil model with additional in- 
formation on bulge sizes, where the ratio between size growth 
via merger and bar instability was self-consistentl y followed for 
each g alaxy in the model (see further details in IShankar et al.1 
1201 laf) . From such a catalog we randomly extracted 100 galax- 
ies with bulge sizes mostly grown via bar instabilities finding 
very similar results compared to our random samples. 

In the right panel we compare model predictions with a col- 
lection of data taken from the Literature, listed in Table 1 and 
discussed in Sect[3] We stress that these data are for active and 
inactive BHs the hosts of which have been recognized as pseu- 
dobulges. These data show that despite the different selections 
and methodologies, they share a similar portion of the Mbh- 
A^buige plane. The most striking feature shown by this collection 
of data is that at fixed bulge mass, most of the objects are mainly 
localized below the classical Mbh-Mb u i ge relation, characterized 
by the grey band, by up to a factor of < 100 and over nearly 
two orders of magnitude in bulge stellar mass. The blue and red 
lines are the predictions for pseudobulges of the B06 and Gil 
models, respectively, with the solid and dashed lines marking 
their respective mean and 1-cr uncertainties (computed from the 
16th and 84th percentile of the logMbh distributions competing 
to each bin of bulge stellar mass). 

We find that the B06 model predicts a rather tight Mbh-Mb u i ge 
relation for all galaxies with normalization consistent with that 
of classical bulges. Thus the model can only marginally account 
for BHs in later type galaxies that observationally mostly lie well 
below the classical Mbh-Mb u i ge relation. The Gil model per- 
forms better in this respect. Although a significant fraction of 
the local BH population clusters around the classical Mbh-Mb u i ge 
relation, a non-negligible portion extends to very low mass BHs 
at fixed Mb u i ge - The solid and dashed blue lines highlight the dis- 
tribution of galaxies in the Gil model that have grown most of 
their bulge mass via disk instabilities. More precisely, we find 
that moving, at fixed Mb u i ge , from the lowest mass BHs to the 
highest ones, the fraction of mass assembled via mergers grad- 
ually increases. The merger model therefore does not predict a 
bimodality of BHs in classical and pseudo bulges, but rather a 
continuous trend. 

A more quantitative comparison between models and data 
is provided in Figure [2] which shows BH mass distributions 
for a fixed bin of bulge stellar mass. Given that we are here 
mainly interested in the comparison among predicted and mea- 
sured scatters in BH mass at fixed bulge mass, we simply nor- 
malize each histogram to have the same peak at 1 . Distributions 
for galaxies with bulge masses 10 < logMb u i ge < 10.5 and 
10.5 < logMb u i g e < 11 are plotted in the left and right panels, 
respectively. The red and blue lines are the predictions for the 
Gil and B06 models, respectively, separated for classical and 



pseudo bulges (dotted and solid lines, respectively). The solid, 
black lines in both panels are the data collected in TableQ] while 
the grey stripes are the b lack hole mass in tervals (median plus 
scatter) inferred from the ISani et alJ d201 lh Mbh-Mb u i ge relation 
for the corresponding mean bulge masses. 

For each log Mbh distribution we then compute the median 
value, and the 16th and 84th percentile which define the ±l<x 
errors on logMbh- We find that, at least above logMb u i ge Z H 
both models predict a narrower distributions for classical bulges, 
with an intrinsic cr(log Mbh) ~ 0.3 - 0.4 dex, in good agree- 
ment with the empirical one of 0.38 dex (grey band). However, 
while the B06 model continues to predict a similarly narrow 
distribution for pseudobulges, the Gil model produces much 
broader distributions for pseudobulges, with an intrinsic scatter 
of cr(log Mbh) ~ 0.6 dex, a factor of about two higher than for the 
B06 model and in better agreement with the data (cfr. solid, red 
lines, the Gil model/pseudo, with the black, solid histograms, 
the data). At lower bulge masses (left panel), both models envi- 
sion similar BH mass distributions and similar differences, with 
the Gl 1 model better matching observations. 

In the Gil model seed BHs (set to 1O 3 M ) grow via mergers 
by orders of magnitude to the level of > 10 5 - 1O 7 M , com- 
parable to the BH masses measured in late-type galaxies. At 
later times, when mergers become rarer, BH growth is blocked 
while host galaxies continue to grow via in-situ star formation 
and bulges grow via disk instabilities. Thus the predicted Mbh- 
Mbui ge relation in the Gil model will inevitably show a rather 
large scatter at fixed bulge mass. 

4.2. Pseudobulges and Classical bulges: Evolution 

In Fig. 3 we show several mean properties for the subsamples 
of pseudobulges of the B06 and Gil models discussed earlier 
(solid black and blues lines, respectively), compared with the 
mean properties of subsamples of classical bulges (red long- 
dashed and dot-dashed lines, respectively). For each object we 
trace back its more massive progenitor in the merger tree and 
record its properties. IShankar et al.l d201 lab found that the sizes 
of pseudobulges are predicted to be a factor of a few lower than 
those of classical bulges at fixed stellar mass. Fig. 2 extends their 
study to other properties. 

4.2.1. B/T versus z 

We find that the bulge-to-total B/T ratios (upper middle) of 
both galaxy populations increase rapidly at very early epochs. 
However, while classical bulges continue to significantly grow 
and turn into ellipticals (B/T~ 1) at late times, the B/T of pseu- 
dobulges quickly flatten out to values B/T< 0.6 after a few Gyrs, 
as a consequence of their significantly different growth patterns 

B 

4.2.2. S FR versus z 

Pseudobulges in both models reside in lower mass haloes at 
z < 2 (lower left), which contributes in explaining the different 
merger histories; it also induces less powerful mergers, thus less 
BH growth at all times. We find in fact that pseudobulges have 
grown > 70% of their final stellar mass via disk instabilities (see 

2 The apparent slight decrease in the B/T for pseudobulges in the 
B06 model is simply an artifact of the averaging, the large intrinsic 
scatter, and the choice of galaxy subset, but not driven by any evident 
physical process. 
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also discussion in B06), while classical bulges of the same mass 
tend to grow mostly via late dry minor mergers (see details in 
IShankar et alj|2011al and references therein). This is also evi- 
dent from their different cosmological star formation rate (SFR; 
lower middle; information on the SFR are only available for the 
Gil catalogue) and specific SFR (SSFR; lower right). The SFR 
in pseudobulges peaks at later times and their SSFR is larger by 
a factor of a few at z < 2.5 with respect to classical bulges that 
become progressively dead with time though still significantly 
increasing in stellar mass (via dry mergers). 

4.2.3. BH mass versus z 

The characteristic BH mass (upper right) in pseudobulges for the 

Gil model (solid blue) is a factor of several lower (> 4 10 

depending on redshift) than those in classical bulges (red dot- 
dashed) since early times, at variance with the bulge masses that 
differ by only a factor < 2 (upper left). The B06 model instead 
predicts that both BH and bulge masses of pseudo (solid black) 
and classical (red long dashed) bulges are similar because they 
always co-evolve. At late times pseudobulges continue forming 
stars in their disks, undergoing disk instabilities, and growing 
their bulges. However BHs in the Gil model suffer severe starv- 
ing due to less numerous and less powerful merger events, thus 
inducing further departure from the classical Mbh-Mb u i ge rela- 
tional We'd like to conclude stressing that what discussed here 
are possibly extreme BHs in pseudo and classical bulges, but the 
models also predict significant fractions of BHs with intermedi- 
ate properties. 

4.3. Late-type bulges: BH demography revisited 

Having compared theoretical and observational relations of BHs 
in classical and pseudo-bulges, we here move on to discussing 
the full demography of the predicted and observed BH popula- 
tions via the BH mass functions. While this type o f compari- 
son was pursued several t imes in the Literature (e.g., Lap i et all 
20061 iMarulli et al.ll2008h . it was never performed before taking 
care of differentiating classical and pseudo bulges. 

To this purpose, we first recompute an empirical estimate of 
the local BH mass function (BHMF). A common way to ob- 
tain an estimate of the BHMF is via convolution of the velocity 
dispersion function (VDF) of local galaxies, with the M^-cr re- 
lation (the convolution is between the VDF and a Gaussian with 
dispersion equal to the intrinsic scatter in the M\,h-cr relation). 
To this purpose, we use th e latest determ i nation of the Hubble 
type-dependent VDF from iBernardi et ail (l2010h . coupled with 
a revised M^-cr relation by Shankar et al. (2011, in prepa- 
ration) obtained from an updated sample of iFerrarese &Fordl 
d2005h . i.e., logM bh = (4.57 ± 0.47) log cr + 8.26 + 0.08 with 
intrinsic scatter r\ = 0.17 + 0.08 (for early-type galaxies), and 
logMbh = (4.99 + 0.38) logo- + 8. 18 + 0.06 with intrinsic scatter 
rj = 0.17 ± 0.06 for later-type galaxies (though similar conclu- 
sions are derived adopting other estimates of the relation). The 
black, dot-dashed line in Fig. 4 is the result of the convolution 



3 We note that each subset of mock galaxies was mostly chosen on 
the basis of being enough massive to host a super-massive BH, and on 
the relative fractional growth via mergers and disc instabilities. We thus 
expect by construction the pseudo and classical bulge subsets within 
each model to yield similar stellar masses at z = (at the factor of 
< 2 level). However, we do not expect the classical and/or pseudobulge 
subsamples between the two model to necessarily yield the same bulge, 
BH, or virial mass at z = 0. 



including all Hubble types, with the dotted lines marking the 
1-cr error bars derived from Monte Carlo simulations that in- 
clude statistical uncertainties on the VDFs and scaling relations 
(including uncertainties in the intrinsic scatter). The estimated 
BH mass function is in go od agreement with previous estimates 
(e.g., IShankar et al.ll2009l grey band) 

However, iFisher & Drorvl d201 ll) recently suggested that 
pseudobulges may be the dominant class of bulges at low 
masses. Moreover, as discussed above, pseudobulges host signif- 
icantly lower mass BHs than expected from the Mbh-o - relation. 
The long-dashed, cyan line in Fig. 4 shows the local BH mass 
function obtained by assuming that all bulges in Sb galaxies are 
actual ly pseudobulges and thus follow the Mbh-o - relation by Hu 
(2008). (We do not claim that th ere is a separate Mbh-o" relation 
for pseudobulges; the iHul (120081) relation is used here just for il- 
lustration.) The solid, cyan line is the extreme estimate in which 
all bulges in Sb and Sa galaxies are pseudobulges. This exer- 
cise shows that accounting for pseudobulges can easily decrease 
the BH mass density by ~ 15%, when correcting only Sb, up to 
~ 35% in the extreme case that all Sb/Sa galaxies are actually 
pseudobulges, while BH number densities decrease by a factor 
of ~ 2 - 3 below < 10 8 M r (these res ults do not depend much on 
the assumed scatter in the Hu (2008) relation). 

The right panel of Fig. 4 compares the local BH mass func- 
tion corrected for pseudobulges (cyan area, obtained by combin- 
ing the cyan long-dashed and solid lines in the left panel), with 
the global BH mass function at z — predicted by the B06 and 
Gil models (long-dashed, black and solid, blues lines, respec- 
tively). It is interesting to note that both models still predict too 
many classical (i.e., more massive at fixed Mb u i g e) BHs in the lo- 
cal Universe, thus tending to overpredict the number of low mass 
BHs, in the range 10 7 and 10 s M , with respect to the empirical 
BH mass function corrected for pseudobulges (below 10 7 M Q , 
the number density continues to be overpredicted by Gl 1, but it 
is underpredicted by B06). 

From this rather straightforward experiment we thus con- 
clude that if BHs in pseudobulges really are the most common 
type of BHs in the local Universe, then models sti ll need to prop- 
erly account for such a large population (see also Fontano t et alj 



2011; Shankar et al 



2011 ah. Recent estimates of the local BH 



mass function (e.g., Vik a et alJl2009T) claim even lower number 
densities of low mass BHs that would further exacerbate the ten- 
sion with model predictions. The shape of the local BH mass 
function at low masses is essential to properly constrain the total 
amount of sub-Eddingt on accretion that cha racterized the cos- 
mic evolution of BHs (Shankar et al. 11201 lbl) . 



5. Discussion 

The comparison in Fig. 1 suggests that a tight link between bar 
instabilities and BH growth at all times and in all conditions may 
be disfavored. More specifically, we have explored the outputs of 
the B06 model that couples BH growth to mergers and disk in- 
stabilities with both fuelling modes tuned to match the median 
local scaling relation between BHs and bulges masses in early- 
type galaxies. We have then compared the model outputs against 
a (recent) collection of data on BH masses in local active and 
inactive galaxies of late-type morphological type, mostly iden- 
tified as pseudobulges, that tend to be outliers with respect to 
the scaling relations of BHs in earlier-type, classical bulges. The 
model struggles to fully reproduce the latter selection of data. 
In other words, the intrinsic scatter predicted by the B06 model 
induced by differences in galaxy-to-galaxy evolutionary paths is 
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Fig. 4. Left panel: estimate of the local black hole mass function derived by convolution of the galaxy velocity dispersion function 
and t he M^h-Q" rela t ion c onsidering all bulges as classical (dot-dashed and dotted lines), compared to the collection of estimates 
from Shankar et al.l d2009l g rey band). The cyan solid line is the estimate of the local black hole mass function assuming that all 
Sbc+Sa galaxies follow the iHul (2009) M bh -o~ relation proper of pseudobulges. The cyan dot-dashed line is the local black hole 
mass function obta ined assuming that on ly Sbc are pseudobulges. Right panel: predicted local black hole mass function from the 
Bow er et"al] ([2006) and iGuo et all d201 lb models (black long-dashed and blue solid lines, respectively) compared with the local 
mass function not corrected (dot-dashed and dotted lines) and corrected for pseudobulges (cyan area, combination of the dot-dashed 
and solid lines in the left panel). 



not sufficient by itself to match the recent data on BHs in later 
type galaxies/pseudobulges. 

The Gil model instead, triggering BHs only in mergers, al- 
lows for significant, late bulge growth via bar instabilities with 
no parallel fuelling in the central BH, thus naturally creating a 
larger scatter in the BH-bulge mass relations. Our results thus 
suggest that models in which late time disk instabilities, mainly 
characterized by the formation of bars, are not closely related to 
BH growth, tend to generate a larger scatter in the Mbh-Mb u i ge 
relation in better agreement with the broad distribution of BHs 
we o bserve today. This finding confirms prior e xpectations (see, 
e.g., FKormend v et al.ll201 U [Mathur et aLl l201 ll and references 
therein). 

We note here that several other galaxy/AGN evolutionary 
mode ls based on BH fuelling in bulges with no bar instabilities 
(e.g. iGranato et al.ll2004t ICirasuolo et"ai1 120051: IShankar et all 
2006) have also suggested additional possible sources of the ob- 
served large scatter observed in BH scaling relations. A break 
in BH scaling relations at low masses induced by inefficiency of 
quasar compared to supernova feedback, r esponsible for shap- 
ing the galaxy stellar mass function (e.g., Sha nkar et al.l 12006, 
and references therein), and different normalizations of the M^- 
Mbuiee relation at different virialization epochs of the host haloes, 
inevitably increase the predicted distribution of BH masses at 
fixed bulge mass. 

Our results, however, do not imply that BH fuelling trig- 
gered by disk instabilities is never at work. Within the highly 
turbulent interstellar medium of the protogalactic disks at high 
redshifts mentioned in Sect. [Tl (e.g.. iGenzel et all 1201 ll) . large 
clumps can form and migrate towards the centre via dynamical 
friction. At variance with bars that can build up bulges gradu- 
ally on very long timescales, the latter type of instabilities are 
rapid and violent processes. Moreover, the clumpy disk insta- 
bilities may potentially form central spheroids structurally more 
similar to early-type bulges, as displayed by some recent numer 



Dlaye 

ical simulations (e.g jBournaud et al.l201 ll but see also Inoue & 



Satoh 201 1). The clumpy inflows within high-z disks have been 
recently recognized to also play a non-neg ligible role in feed- 
ing the central BHs dBournaud et al.ll201 1|). possibly under th e 
effect of continuous cold flow accretion (Di Matt eo et al.ll201 ll) . 
In this framework, the turbulent disks at high redshifts seem to 
share similar effects to mergers, formin g classical bulges and 
proportionally growing central BHs (e.g. jBournaud et aDl201 ll 
and references therein). 

Our results may equivalently reflect on a still not fully ap- 
propriate treatment of BH fuelling during bar-like disk instabil- 
ity events. The way BH feeding follows late bar instabilities may 
in fact be a non-trivial function of other parameters such as gas 
fractions, environment, etc... that may further contribute to pro- 
duce the local scatter we observe. 

Our understanding of BH feeding is in fact still fuzzy, espe- 
cially at lower redshifts when gas reservoirs in galaxies become 
poorer. Peeples & Martini (2006) raised some doubts about a 
strong link between bars and AGN activity. They investigated the 
correlation between bar strength and circumnuclear dust mor- 
phology for 75 galaxies finding that in many cases strongly 
barred galaxies showed structures terminating in a circumnu- 
clear starburst ring at ~ 10 pc away from the galactic nucleus. 
On the other hand, many Narrow Line Seyfert 1 galaxies appear 
as barred systems and their AGNs are believed to be fuelle d via 
bars (e.g JCrenshaw et al.ll2003t lOrban de Xivrv et al.ll201 ll) . 

ICoelho & Gadottil (1201 ll) recently claimed evidence for dou- 
ble the fraction of AGNs in barred galaxies as compared to un- 
barred galaxies, at least for low-mass bulges. lOh et all (1201 ll) 
have recently analyzed a large sample of barred and unbarred 
galaxies extracted from the Sloan Digital Sky Survey conclud- 
ing that the correlation between bars and AGN activity (and 
SFR) is far from trivial and also difficult to establish due to sec- 
ondary correlations between AGNs, galaxy morphology, etc... 
They however point out that some connection between AGN ac- 
tivity and bars seems to be present, at least in the local Universe. 
In particular, they claim that galaxies having longer bars appear 
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to have more AGNs than their counterparts at fixed stellar mass 
and morphology. 

Indirect evidence for BH growth via ba rs may also come 
from higher redshifts. ICisternas et al.l (1201 ll) find the BH to to- 
tal host stellar mass ratio to be constant since redshift z = 0.9. 
Given that their sample is mainly composed by galaxies with 
disks, they claim that secular processes might have fuelled 
both a central bulge and correlated BH in time thus preserv- 
ing their correlation down to z = 0. More generally, other ad- 
ditional sources of BH fuelling at later time s are also proba- 
bly a t work, such as winds from stars (e.g., Ciotti & Ostrikerl 
2001; Kauffma nn & Heck man 200l ICenll20T H and references 
therein). 

Overall, evidence for some AGN activity in 
barred/pseudobulged systems may be present. However, 
there is also mounting evidence that these systems lie mostly 
below the relations characteri zing classical bulges , for both 
active and inactive galaxies. iBeifiori et al.l (1201 2|) have re- 
cently revisited a number of BH scaling relations in a revised 
sample of local galaxie s with BHs, and claimed, in line with 
Kormend v et al.l (1201 ll) . that pseudobulges have little or no 
correlation between BH mass and velocity dispersion. Thus, it is 
not at all certain that BHs in n on-classical bulges sh o uld en d up 
obeying similar scaling laws. lOrban de Xivrv etall (1201 ll) . for 
example, suggest that Narrow Line Seyfert 1 may not be in any 
special phase of their evolution, but simply are BHs growing 
slowly due to their specific duty cycle and spin, mainly driven 
by secular processes in their ps eudobulge hosts. 

As discussed in Sect. 14.31 Fisher & Drorvl (1201 ll) recently 
suggested that pseudobulges may be the dominant class of 
bulges at low masses. Thus, the number of outliers in the BH- 
galaxy classical relations could be very high, impacting the de- 
mography of local BHs by severely reducing the expected num- 
ber density of low mass BHs in the local Universe (similar 
thoug hts were also pu t forward by Gadotti & Kauffmann 2009 
and lGreene et al.l2010t) . This in turn could have a non-negligible 
repercussion on accretion and merging models that try to tune 
parameters such as the radiative efficiency and Eddington ratio 
distributions b y matching the full shape of the local BH mass 
functio n (e.g.. lShankarl2009HMarulli et al.ll2008l;IShankar et al. I 
1201 lbl) . 



Table 1. Galaxies hosting "pseudobulges" with known black 



hole masses 


Galaxy 


Log M bh [M ] 


Log M buhx [M Q ] 


Reference 


NGC1068 


6.92 


10.60 


H09 


NGC3079 


6.40 


10.24 


H09 


NGC3393 


7.49 


10.62 


H09 


Circinus 


6.04 


9.62 


H09 


IC2560 


6.46 


10.31 


H09 


MW 


6.61 


10.11 


H09 


MW 


6.63 


10.28 


Kll 


Circinus 


6.08 


10.14 


Kll 


NGC1068 


6.93 


11.20 


Kll 


NGC1300 


7.85 


10.14 


Kll 


NGC2748 


7.67 


9.80 


Kll 


NGC2787 


7.61 


9.80 


Kll 


NGC3227 


7.18 


9.97 


Kll 


NGC3384 


7.26 


10.52 


Kll 


NGC4736 


6.82 


10.42 


Kll 


NGC4826 


6.13 


10.30 


Kll 


NGC7582 


7.74 


10.29 


Kll 


MS2254 


6.60 


10.25 


Mil 


RXJ1209 


6.75 


10.53 


Mil 


RXJ1117 


7.33 


10.49 


Mil 


RXJ2217 


7.10 


10.45 


Mil 


RXJ1702 


7.34 


10.53 


Mil 


KUG1136 


6.89 


9.91 


Oil 


MRK0042 


6.22 


9.91 


Oil 


Mrk0359 


6.70 


10.65 


Oil 


Mrk0382 


6.59 


10.26 


Oil 


Mrk0493 


6.35 


10.07 


Oil 


Mrk0766 


7.14 


9.61 


Oil 


Mrk0896 


7.01 


10.25 


Oil 


Mrkl044 


7.20 


10.03 


Oil 


NGC4748 


6.95 


10.07 


Oil 


Mrk0871 


7.67 


10.22 


Oil 


Mrkll26 


6.90 


9.80 


Oil 


Circinus 


6.23 


10.26 


Sll 


IC2560 


6.64 


10.94 


Sll 


NGC1068 


6.93 


11.33 


Sll 


NGC3079 


6.40 


11.02 


Sll 


NGC3368 


6.88 


10.74 


Sll 


NGC3489 


6.78 


10.07 


Sll 


NGC3998 


8.38 


10.36 


Sll 


NGC4258 


7.58 


10.98 


Sll 


NGC4594 


8.76 


11.23 


Sll 
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Notes. Details on the data listed in this Table can be f o und i n 
Sect. [3] Here H09, K ll Mil, Oil a nd SI 1 refer to iHul (2 009), 
iKormendv et alJ <201lJ). iMathur et alj d201ll) . lOrban de Xivrv et all 
j201 ll) . and lSani et al.1 ( 1201 ll) . respectively. 
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